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Charm- and bottom-quark rescattcring in a Quark-Gluon Plasma (QGP) is investigated 
with the objective of assessing the approach toward thermalization. Employing a Fokker- 
Planck equation to approximate the collision integral of the Boltzmann equation we augment 
earlier studies based on perturbative parton cross sections by introducing resonant heavy- 
] light quark interactions. The latter are motivated by recent QCD lattice calculations which 

■ indicate the presence of "hadronic" states in the QGP. We model these states by color- 
' less (pseudo-) scalar and (axial-) vector D- and _B-mesons within a heavy-quark effective 

(N . theory framework. We find that the presence of these resonances at moderate QGP tcmper- 

^ I atures substantially accelerates the kinetic equilibration of c-quarks as compared to using 

perturbative interactions. We also comment on consequences for D-meson observables in 
ultra-rclativistic heavy-ion collisions. 

PACS numbers: 12.38.Mh,24.85.-hp,25.75.Nq 
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> ' I. INTRODUCTION 

in 

■ Hadrons containing heavy quarks are valuable probes of the strongly interacting matter pro- 
^ , duced in high-energy collisions of heavy nuclei. The spectral properties of (bound) c-c states, 
^ I such as their binding energy and (decay) width, are expected to undergo substantial modifications 

' in a Quark-Gluon Plasma (QGP) and thus affect charmonium yields and momentum spectra in 

heavy-ion reactions with QGP formation (see, e.g., Refs. Q, y] for overviews). Reinteractions 
of individual c-quarks in the QGP will reflect themselves in transverse-momentum (pt-) spectra 
of open charm hadrons (L'-mesons) PiSSQ; ™ost notably their elliptic flow, V2{pt), in semi- 
^ . central collisions 0, 0] . Preliminary experimental results from the Relativistic Heavy-fon Collider 

^ ' (RHIC) indicate the possibility that the Z)-meson V2 could be similar in magnitude to the one of 

light hadrons 3, Since the c-quark is rather heavy, this would be quite remarkable and could 



^ . provide important insight into (nonperturbative) properties of the QGP at moderate temperatures, 

^1 T ~ 1-2 Tc- E.g., in the light-quark sector, parton rescattering through hadron-like states in the 

QGP (motivated by recent lattice calculations of Quantum Chromodynamics (QCD) at finite tem- 
perature) has been suggested as a mechanism to enhance partonic cross sections ESQ in order 
to facilitate rapid thermalization of the bulk matter at RHIC as required in hydrodynamical mo- 
dels. The notion of charmonium resonances in the QGP Hill has been applied earlier to assess 



J ftp production at SPS and RHIC Employing kinetic rate equations to consistently account 
for both dissociation and regeneration reactions, c + c<^ J/tp + X, J/ip resonance formation in the 
QGP via c-c "coalescence" turns out to be the dominant contribution to the final yield in central 
Au-Au collisions at RHIC 12, 3,^^ (see also Refs. |2fll. l2]]|l. Furthermore, J/V' production was 



found to be sensitive to the in-mediuni properties of charm quarks, i.e., their in-medium masses 
and degree of kinetic equilibration ESQ. 

Early studies of charm-quark thermalization in the QGP have been conducted in Ref. based 
on elastic perturbative QCD (pQCD) cross sections |2j], c + q{q,g) — > c + q{q, g), implemented 
into a Fokker-Planck equation to approximate the collision integral of the Boltzmann equation (see 
Ref. for a recent application to light partons). With a strong coupling constant as=0.6 rather 
short relaxation times of around ~4 fm/c have been found for a massless QGP at temperatures 
T~400 MeV. However, since the relaxation times are essentially proportional to a""^ , more mod- 
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erate values of Os {e.g., 0.3) lead to a significant increase (factor 3-4), rendering thermalization 



Q 



of c-quarks under RHIC conditions unlikely |27l |. This is also consistent with recent transport 
studies Sii]. 

In the present article we evaluate heavy-quark rescattering in the QGP via "D"- and "i?"- 
meson resonances, the existence of which is the main assumption in our work. Although lattice 
QCD has not yet addressed in-medium heavy-light (Q-q) spectral functions at finite temperature, 
resonance-like correlations in the QGP are quite plausible in view of the indications for both q- 
"iol . iH and Q-Q systems. Further support for this assumption is pro vided by calculations using 
effective four-quark interactions within Nambu- Jona-Lasinio models |32l . IH^ . |3] • With form factors 
and coupling constants adjusted to free D-meson masses, finite-temperature calculations lead to 
resonances above the Q-q threshold in the QGP, with appreciable widths of several hundred MeV. 

Our article is organized as follows: in Sec. ^ we introduce Lagrangians for resonant c-q inter- 
actions based on chiral and heavy-quark symmetry for pseudo- / scalar and axial- /vector multiplets 
fSec. lII a)) and evaluate pertinent scattering amplitudes (Sec. Ill B)) . In Sec. lIIII we apply the result- 
ing cross sections within a Fokker-Planck equation; we first determine temperature and momentum 
dependencies of drag and diffusion constants in a static QGP (Sec. IIII A|) . and then evaluate the 
time evolution of c-quark transverse-momentum (pt) spectra within an expanding fireball model to 
simulate conditions in heavy-ion collisions at RHIC (Sec. nil's]) . We conclude in Sec. IIVI including 
a discussion on open-charm observables in ultrarelativistic heavy-ion collisions (URHICs). 

II. D-MESONS IN THE QUARK-GLUON PLASMA 
A. Heavy-Light Quark Lagrangians 

Our description of D-meson resonances in the QGP is based on a rather simplistic quark-meson 
model, accounting, however, for the relevant symmetries, i.e., chiral symmetry in the light quark 
sector (u-d) and heavy-quark symmetry for c-quarks. As a minimal set of resonances, consistent 
with lattice calculations we assume the lowest-lying pseudoscalar (D) and vector mesons (D*) 
to survive above Tc^. In addition, (approximate) restoration of chiral symmetry mandates the 
existence of pertinent SU{2)f chiral partners in the scalar (-Dg) and axial-vector (-Di) channel with 
mass and width identical to D and D*, respectively. The effective Lagrangian thus takes the form 



(1) 



^Dc, + ^l, - iGs (q^^c - qj'^^^c + h.c. 

- Gv (g-T'^^i^^c - q^^j^^.^^c + h.c)j 

with the usual free (kinetic and mass) terms for quarks and D-mesons, 

= c(i^ - mc)c + qi^q, 

The fields $ represent anti- £)-mesons, transforming as isospinors under isospin rotations. 

The interaction terms in Eq. (fTl) will be evaluated to leading order in 1/mc according to heavy- 
quark effective theory (HQET) |35|. This ensures the absence of unphysical (4-D longitudinal) 



In the vacuum these are associated with _D+(1870), I>"(1865) and I>*(2010) mesons. 
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degrees of freedom for massive (axial-) vector meson fields as encoded in the transversality con- 
straints 

v^^^f' = v^^'t = (3) 

{Vfj^: four velocity of the charm quark or D-meson). Equivalent relations hold for the D-meson 
self-energies (to leading order in the l/rric expansion of HQET). 

In the light-flavor sector the relevant symmetry is the invariance of the Lagrangian under the 
chiral SU{2)l x SU{2)ji group. For infinitesimal angles S4iv,A, it is characterized by the standard 
vector and axial-vector transformations acting on the quark fields as 

g ^ (1 + i6(j)vt + iS(j)Atj5)q, c (4) 

with t = t/2 the generators of SU{2) (r: Pauli matrices). To construct the group operations for 
the D fields, we identify the latter with underlying quark currents according to 

$ ~ C75g , ^>o ~ eg , ~ cjf,q , ~ n^j^q . (5) 

Here, "~" denotes "transforms under (jlj like". Thus the transformation rules for the l)-meson 
fields follow as 

$ ^ $ + i6(j)A ■ i"^ + ^S(j)v ■ t^O' 
^>o ^ $0 + i'^'^A • + i5(t)v ■ t<^ 

With these properties the Lagrangian, Eq. (^3), is a scalar under chiral transformations. 

For D^-mesons we restrict ourselves to the (experimentally known) pseudoscalar and vector 
states since spontaneous chiral symmetry breaking in the strange-quark sector is expected to per- 
sist to temperatures significantly larger than Tc (characterized by sizable values of the strange 
condensate {ss)). The corresponding Lagrangian is therefore taken to be 

- iG,,5^-7'^s^c - Gsys^'^l^^^c . 

Eqs. and ((7j) constitute our basic vertices for effective charm-light quark interactions in the 
QGP (analogous expressions hold in the 6-quark sector upon the replacement c — > 6 and D ^ B). 
The underlying parameters, i.e., the bare masses of the D-meson and their coupling strength to 
the quarks, will be fixed to resemble more microscopic model calculations of corresponding spectral 
functions above Tc, as discussed below. In the following we will also impose the spin symmetry of 
HQET implying that the spectral properties of pseudoscalar and vector mesons, and consequently 
their masses and couplings to quarks, are equal. 



(6) 



(7) 



B. Meson Self-Energies and Heavy-Light Quark Scattering Amplitudes 

The key ingredient for the heavy-quark scattering amplitudes are the heavy-meson exchange 
propagators, 

which are essentially determined by the underlying self-energies, Hd b-, together with the bare 
resonance masses, tud b, of the Lagrangian. In the following, we will evaluate the self-energies 
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in terms of the heavy- light-quark loop (cf. the left diagram in Fig. thereby investigating two 
different schemes to regularize the divergent loop integrals to assess the robustness of our results. 
We will concentrate again on the charm-quark case, but completely analogous expressions apply 
to the bottom sector. 

Within the dimensional regularization scheme (see, e.g., [3)) the interaction vertices of Eq. 
yield a self-energy for (pseudo-) scalar L)-mesons of four-momentum k of the form 



3G2 Ami 



2s 



87r2 



4-, 

4 

H ^In 

s 



+ 3ml -2s + {s- 2m, 



7 + In 



m% 



47r^f^ 



mi 



(9) 



(s = A;2). The first term contains the quadratic divergence for d — > 4, showing that the self-energy 
can be rendered finite with a field- and mass-renormalization; fi denotes the mass scale in the 
dimensional-regularization scheme to keep the momentum dimensions of the integrals as for d = 4, 
and 7 ~ 0.577 is Euler's constant. The (regularization-independent) imaginary part is given by 



ImnD(s) 



3G| {s-mlf 



@(s-mj.) . 



Sir s 

For the vector and axial- vector mesons, we employ the HQET propagator for the charm quark 

rric 1 + ^ 



(10) 



mcV ■ I + ir] 2 



(11) 



to obtain a transverse self-energy in leading order of the expansion in l/rric. In Eq. the 
"residual momentum" I of the charm quark is defined via its total four- momentum as Ic = mcV + 1. 
For the D-meson fields the corresponding decomposition for the heavy-quark expansion reads 



kc = mcV + k, s = k^ 



m^ 



1 + 



2v ■ k 



rrir 



V ■ k 



2mr 



+ (12) 



mt 



Applying the Feynman rules, and after some (Dirac) algebra, the dimensionally regularized pola- 
rization tensor for the axial-/vector becomes 



2 



mr 



{2-kY [{I +pY + ir/] {mcV -l + ir]) 
The integral is conveniently evaluated with help of the identity js^] 

2 

(a + 2A6)2 ■ 

Upon integrating over A one obtains for the imaginary part 



1 

ab 



dX- 



imn^':(p) = imn^';(p) 



/ N^Gy (s — m^)^ „ , 2\ 
-{vij.Vu - 9tiu)-^ — B(s - mj . 



8tt 



mf. 



Since up to corrections 0{vk/mc) one can identify mf. = s in the denominator, one finds 



(13) 



(14) 



(15) 



(16) 



5 



which is the expected result from spin symmetry of the HQET. We define the renormahzation 
constants for the self-energy by the following conditions, 

a.n£^")(s)|,=o = o, Reng^'^)(.)|,=^.^=0. (17) 

The first condition ensures that, within a vector dominance model, the photon propagator has 
residuum of unity at s = 0, while the second one implies that the renormalized meson mass coincides 
with the bare mass of the Lagrangian. The renormalized self-energy is, of course, independent of 
the dimensional-regularization scale fx. 

As an alternative way of regularizing the divergent self-energy integrals we introduce a dipole 
form factor at the c-q-D vertex, 




where q denotes the three-momentum of the quarks in the center-of-mass frame. The imaginary 
part of the self-energy is then given by 

Imn^^\s) = lmUD{s)FWq\) , (19) 

\q\ = {s — ml) / {2y/s), while the real part is determined by an unsubtracted dispersion integral. 
The bare meson mass is then adjusted to render a vanishing real part of the propagator at the 
physical resonance mass. 

As default parameters for our calculations we use massless light quarks, a charm-quark mass 
of rric = 1.5 GeV and physical D-meson masses {rnF^^^)"^ = m'j^ — Re n£)[(m^-^'')^] = 2 GeV 
(corresponding to a vanishing real part in the propagator). The coupling constant G = Gsy is 
varied to allow for widths of the D-meson spectral functions of 300-500 MeV, to approximately 
cover the range suggested by effective quark models • It is important to note that we 

assume the D-meson resonances to be located above the c-q mass threshold, mc+rrig, which renders 
them accessible in c-q scattering processes. The situation is quite different for (bound) meson states 
(i.e., below the anti-/quark threshold), where the resonant part of the scattering amplitude cannot 
be probed through c + q — > c + q interactions (even for resonance masses close to threshold, thermal 
energies of anti- / quarks imply that the average collision energy is significantly above the resonance 
peak). In this case, other processes need to be calculated, e.g., c -|- g — > D + g, where the extra 
gluon in the final state carries away four-momentum to allow the D-meson to emerge on-shell. 
The same framework is also applied to the bottom sector, with 6-quark and i?-meson masses of 
nih = 4.5 GeV and tub = 5 GeV, respectively. 

As we will see below, for equal on-shell masses and widths of the L'-meson spectral functions, 
both regularization schemes lead to quite comparable results for the thermal relaxation properties 
of c-quarks, even though the off-shell properties differ significantly. 

The D-meson propagators figure into the invariant matrix elements for elastic c-quark scattering 
off quarks, c + q ^ c+ q, and antiquarks, c + g — > c + g, in the u- and s-channel, respectively. One 
finds 

l-^gf = 720G\s - mlf |Dd(s)|' (20) 
\M,\'' = 720G^(n - mlf \DD{ut , (21) 

where we have summed over the contributions of all light-quark resonances, being equal due to 
either the heavy-quark symmetry (D-D*, Dq-Di) or chiral symmetry {D-Dq, D*-Di). We also 
included finite-mass strange quarks along with Dg- and D*-mesons. 
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FIG. 1: Left panel: c{b)-q loop diagram representing the £'(i?)-nicson self-energy in the QGP. Right panel: 
" meson" -exchange diagrams contributing to the invariant matrix elements for the scattering of charm quarks 
on light quarks (u-channel) and anti-quarks (s-channel). 



In addition to resonant interactions, elastic scattering in pQCD is accounted for to leading order 
in as, 0{a1). The corresponding matrix elements [2^ have been supplemented with an additional 
gluon Debye (screening) mass, fig = gT, which regulates the forward singularity in the t-channel 
exchange graphs E^?]. The strong coupling constant will be varied over the range ag =0.3-0.5. 



III. CHARM QUARK RESCATTERING IN THE QGP 
A. Fokker-Planck Equation, Drag and Diffusion Coefficients 

The above matrix elements are now implemented within a kinetic theory framework to assess 
the thermalization time scales for heavy quarks in a QGP. Following the steps outlined in Ref. [23 ]. 
we start from the Boltzmann equation for the heavy-quark distribution function f{t,p) and neglect 
any mean-field terms. Furthermore assuming the scattering processes to be dominated by small 
momentum transfers one arrives at a Fokker-Planck equation describing the time evolution of / in 
momentum space, 



df{t,p) ^ d 
dt dpi 



fit,p) . 



(22) 



For an isotropic (rotationally invariant) plasma the drag and diffusion coefficients in (|22() can be 
decomposed as 



Mp) = piAi\p\) 



with the scalar functions 



Bijip) = 

Aim = (1) 

Bo{\p\) = 

Bim = 



p2 



{p ■ p' 
p2 



t/2\ 



Bo{\p\) + ^B,i\p\) 



{{p-p'f 
p^ 



{{P-P'f 



p2 



2{p'-p)+f{l) 



(23) 
(24) 

(25) 
(26) 

(27) 
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FIG. 2: Upper panel: drag coefRcient A (left) and diffusion coefficient Bq (right) as a function of c-quark 
three-momentum at a temperature of T = 200 MeV for various values of as (pQCD scattering) and Ja- 
meson widths (resonance exchanges). Lower panel: The same quantities as a function of temperature at 
fixed three-momentum \p\ = 0. 



The averaging is defined by 

X {2TTf6^^\p + q-p' -q')f{q)X{p') , (28) 

where p {p') and q {q') denote the momenta of the incoming (outgoing) charm- and light- 
quark/gluon, respectively, and f{q) are the thermal Maxwell-Boltzmann distribution functions 
of the light partons^. After integrating over the four- momentum conserving (^-function, and mak- 
ing use of Lorentz invariance of the matrix elements, the expressions for the scalar coefficients, 



^ As pointed out in Ref. [23 the use of quantum distribution functions (Bose-Einstein and Fermi-Dirac) induces 
moderate corrections to the drag and diffusion coefficients obtained from pQCD rescattering; this is still true for 
the "D" -meson resonance interactions as they involve contributions from thermal anti-/quarks only which, at the 
temperatures under consideration, are mostly non-degenerate. E.g., at r=300 MeV, the deviations when using 
Fermi distributions amount to ~20%, which is inside the uncertainties associated with the resonance parameters 
as discussed below. For the total coefficients (including both pQCD and resonance rescattering) the effects of 
quantum statistics approximately cancel. In the following, we will use Maxwell-Boltzmann distributions for light 
partons, since (i) it is more consistent within the Fokker-Planck treatment for the heavy-quark, the thermal limit 
of which is a Boltzmann distribution, and (ii) it minimizes the deviations in the dissipation-fluctuation theorem, 
cf. the discussion following Eq. 14511 below. 
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FIG. 3: Left panel: Comparison of the friction coefficient for resonances with self-energies calculated in the 
renormalization scheme ()17|l (solid line) and with form- factor regularization (dashed line). The coupling 
was chosen such that in both the resonance width r=0.4 GeV at the resonance mass. Right panel: The 
friction coefficient for different resonance masses as function of the temperature, evaluated in the form-factor 
regularization scheme. 



Eqs. H25|) - (|27|) . can be reduced to numerically tractable three-dimensional integrals'^. The temper- 
ature and momentum dependencies of the A and Bq coefficients are summarized in Fig. |^ The 
main finding here is that for resonant rescattering both are increased by a substantial factor (~3) 
over the pQCD results (the latter are very similar to Ref. [3] as the extra factor of 2 is essen- 
tially compensated by the larger screening mass, fig = gT, used in our calculation). The main 
reason for this effect is not so much an increase in the total cross section (the peak cross section 
in the resonance case of about 10 mb is comparable to an almost constant 4 mb for pQCD), but 
the isotropic angular distribution of the resonance cross sections in contrast to forward dominated 
pQCD rescattering. Also note that the coefficients are rather insensitive to the underlying coupling 
constants for both resonance and pQCD scattering; naively, the matrix elements are proportional 
to the 4th power of the coupling constant, which would imply a variation by a factor of (5/3)^ ~2. 6 
for the coefficients A and Bq for the parameter ranges shown in Fig. [21 The much smaller actual 
variation is due to compensating effects induced by an increased Debye mass for pQCD t-channel 
gluon exchange, and by an increased resonance width for s-channel D-meson exchange. 

To further illustrate the uncertainties inherent to the resonance properties, we display in Fig. |31 
the sensitivity of the drag coefficient with respect to the regularization schemes and the resonance 
masses (for form- factor regularization). From the left panel one observes that the momentum 
dependence of the A coefficient is somewhat more pronounced for the form-factor regularization, 
but the absolute magnitude in both schemes is very similar. The right panel indicates a more 
pronounced sensitivity to the value of the resonance mass. In particular, it confirms that D-meson 
masses close to threshold lead to a significant reduction of the drag effect, due to the thermal motion 
of the light partons from the heat bath. For the same reason, larger masses imply a steeper increase 
of A with temperature (of course, above T ~ 2Tc the very existence of resonance correlations is 
questionable) . 



^ We note that, wlien using tlie pQCD matrix elements, our final expressions are larger by a factor of 2 than the 
ones originally derived in Ref. 1241 . which confirms the findings of Ref. |27|. 
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B. Time Evolution of Momentum Spectra 

1. Time Dependence of the Fokker Planck Equation 

To obtain a better estimate of the effect of heavy-quark rescattering on the transverse- 
momentum (pt) spectra of heavy-flavor hadrons in URHICs we investigate in this section the 
time evolution of the Fokker-Planck equation in a thermally evolving QGP. The latter is modelled 
by an expanding fireball under conditions resembling central Au-Au collisions at RHIC. The main 
simplifying assumption consists of momentum-independent drag and diffusion coefficients. Ac- 
cording to the discussion at the end of the previous section, we will therefore employ D-meson 
self-energies within the renormalization scheme H17() (with masses and widths of m£)=2 GeV and 
r=0.4 GeV), as well as pQCD cross sections (with a<,=0.4), cf. Fig. El The Fokker-Planck Eq. ^ 
then takes the form 



df d d' 



l{t)^{m+D{t)^f , (29) 



where 7 = A{T{t), \p\ = 0) and D = Bo{T{t), |^ = 0) = Bi{T{t), \p\ =0). The time dependence 
of the coefficients enters through their dependence on temperature (as determined in the previous 
section), with T[t) following from the fireball model outlined below. For the initial condition, 

/(t = 0,p)^/o(p1 , (30) 

we will employ c-quark spectra from p-p collisions as extracted from the PYTHIA event gener- 
ator [i^ . lig! ]. The initial- value problem can be conveniently solved employing Green's function 
techniques: if we can find a solution G(t,p;po) to Eq. (|29|) with the initial condition 

G{t = Q,p-po) = 5^''\p-po) , (31) 
the full solution with an arbitrary initial condition l|3()j) follows as 

f{t,p) = j d3poG(t,p;#o)/o(#o) . (32) 

To determine the Green's function, we define its Fourier transform, 

Git,p;po) = j Sq exp{-iq-p) git,q;po) , (33) 

and insert it into Eq. ()29() . leading to the first-order differential equation 

I- -I (34) 

With the initial condition for g determined by Eq. H31|l . 



5(0, q,Po) = j^-:^^ exp(-ipo • q) , (35) 



its solution reads 



g{t,q;Po) = exp {ipo ■ qexp[-T{t)]} exp[-A{t)q^] , (36) 
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where 



dT7(T) 



A(t) = exp[-2r(t)] [ dTZ)(r)exp[2r(T)] . 
Jo 

The Fourier transformation (|33j) yields the result for the Green's function, 

(p-poexp[-r(t)])2 



Git,p;po) 



1 


3/2 . 


exp <^ 


_47rA(t)_ 



4A(t) 



(37) 
(38) 

(39) 



and the time evolution of the distribution function, Eq. I|32|) . is finally obtained by numerical 
integration. 



2. Limiting Cases 

Before we turn to the full solution, let us first illustrate a few limiting cases. For time- 
independent coefficients 7 and D, Eqs. (jS7)l and (P5]l simplify to 



r(t) = 7t , A(t) = ^ [1 - exp(-270] . 
27 



(40) 



Inserting these into Eq. ()39() yields 
Git,p;po) -- 



7 1^^^ ^ f 7 [P -Poexp(-7t)]^ 

27rL»[l - exp(-27t)] J \ 2D 1 - exp{-2'yt) 



(41) 

which was already derived in Ref. |2^. This, in particular, shows that 7 has the meaning of a drag 
(or friction) coefficient, 

(p(t)) =Poexp(-7t) , (42) 
characterizing the equilibration time scale r = I/7. With the momentum fluctuation given by 

3D. 



{p^it))-my = —[i-eM-27t)] 

7 



(43) 



D is readily identified as a momentum diffusion coefficient. From the left panel of Fig. |3]we see 
that for temperatures expected in central Au-Au collisions at RHIC, the equilibration time scale for 
charm quarks is substantially reduced by resonant interactions, to about a few fm/c, comparable 
to the duration of the (putative) QGP phase. Even though a similar mechanism is operative 
for bottom quarks, their much larger rest mass renders the pertinent kinetic equilibration time 
significantly larger, around 10 fm/c or more. 

For t ^ 00, Eq. (|4H) approaches a Maxwell-Boltzmann distribution, 



lim G{t,p;po) = fcq{p) 



27rD 



\3/2 


r —2 n 

IP 


j exp 


2D 



(44) 



and thermal equilibrium implies the dissipation-fluctuation theorem. 
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FIG. 4: Left panel: equilibration time-scale t = I/7 for charm and bottom quarks in the QGP as function of 
temperature with (solid and dashed line) and without (dash-dotted and dash-double-dotted lines) in-medium 
resonances. Right panel: consistency check of the dissipation- fluctuation relation, Eq. H45|l. for c-quarks with 
(solid line) and without (dash-dotted line) resonances, as well as for 6-quarks (with resonances, dashed line), 
in the QGP. 



The such obtained temperature T should, of course, coincide with the one of the (light-quark and 
gluon) heat bath entering through the thermal distribution functions in Eq. H28|) . Thus, Eq. (jlSJ 
serves as a consistency check for the determination of the drag and diffusion coefficients within our 
model, especially for the assumption on the dominance of small momentum transfers underlying 
the Fokker-Planck Equation (|22|) . From the right panel in Fig. |3]we see that, for charm quarks, 
the dissipation-fluctuation theorem is well satisfied (within 3%) when using forward peaked pQCD 
cross sections; but also for the isotropic resonance cross sections within the renormalization scheme 
the deviations do not exceed 11% even at the highest considered temperatures, while within the 
form-factor regularization scheme they reach up to 26% (the latter value is reduced to ~17% when 
an average over a thermal momentum distribution is performed). The latter is due to the greater 
variation of the friction and diffusion coefficients with momentum which makes the approximation 
of momentum-independent coefficients, underlying the derivation of the fluctuation-dissipation 
relation ()45|). less accurate. For this reason, in the next section we shall use the renormalization 
scheme without form factor to investigate the time evolution of pT spectra. For the heavier bottom 
quarks the fluctuation-dissipation theorem is satisfied to high accuracy. 



3. Charm-Quark pT-Spectra at RHIC 

Let us finally address the time evolution of charm-quark pj^-spectra including the temperature 
dependence of drag and diffusion coefficients. To obtain the time evolution of the temperature we 
use a simple expanding fireball model [i^]. In reminiscence to hydrodynamic simulations |4ll | of 
central Au-Au collisions at RHIC, the fireball volume is parameterized by 

VfbW = 7r(zo + v,t){rQ + ]^a^t^f, (46) 

where rg = 6.5 fm and zq = 0.6 fm are the initial transverse and longitudinal size (the latter 
corresponding to a formation time of tq = 0.33 fm/c). The longitudinal and transverse expansion 
are characterized by Vz = lAc (covering a thermal width of about 1.8 units in rapidity) and 
a_L = 0.055 c^/fm (yielding a total fireball lifetime of about 14 fm/c with a thermal freezeout 
temperature of ~110 MeV). 
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FIG. 5: Left panel: results for the time-evolved c-quark pr-spectra in the local rest frame with a temperature 
profile corresponding to QGP and mixed phase in central Au-Au collisions at ^sat jv=200 GeV; dashed curve: 
initial spectrum taken fromp-p collisions; dash-dotted curve: final spectrum using pQCD cross sections only; 
solid curve: final spectrum using both pQCD and D-meson resonance interactions. Right panel: explicit 
time evolution in time steps of 1 fm/c for the pQCD-|-resonancc interactions. 



Assuming isentropic expansion, the temperature at each instant is calculated from the total 
entropy of produced particles, S = s(r)T^B(i) — 10^ (within Ay=1.8), with the entropy density 
in the QGP given by 

s = —T'il6 + W.5Nj) , (47) 

{Nf is the effective number of quark flavors, taken to be 2.5). 

With these parameters, the initial QGP temperature is Tq ~ 375 MeV, decreasing to the 
critical temperature of Tc — 180 MeV after about 3 fm/c, with further evolution in a hadron- 
QGP mixed phase for another 3 fm/c. For simplicity, we treat the latter phase as a QGP at 
constant T = Tc ~ 180 MeV, but accounting for lower parton densities as estimated from the 
temperature dependence of the Fokker-Planck coefficients (from the lower panel of Fig. [21 we 
estimated A,B(x g^^^, where g denotes the density of the light partons). 

For the initial distribution foip), Eq. (|30|) . we employ charm-quark p^-spectra as generated in 
proton-proton {p-p) collisions at 200 GeV by PYTHIA |3g] . A suitable parameterization thereof is 
given by |0 

dpi (1 + pt/BY 

with A = 0.5 GeV, B = 6.8 GeV, a=21, and C=0.845 GeV'^. 

As detailed above, the time evolution of the pT^-distribution is obtained from Eq. ()32() with 
the Green's function (|41|) . integrated over pz, which amounts to a two-dimensional solution of the 
Fokker-Planck equation. Fig. [S] shows the time evolution of c-quark p^^-spectra in central Au- 
Au at RHIC, evaluated in the local (thermal) rest frame of the expanding matter in QGP and 
mixed phase (i.e., the additional boost from the collective transverse expansion is not included), 
using the drag and diffusion coefficients computed in Sec. IIII Al From the left panel one finds 
that, when allowing for pQCD rescattering only, the initial spectra from p-p collisions are affected 
rather little. On the contrary, when augmenting the interactions with Z)-meson resonances^, the 



* Note that, when adding the matrix elements for pQCD and resonance interactions, we did not resum the pertur- 
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pr-spectra undergo a marked reshaping, essentially a redistribution from high to low pT- The final 
spectra (i.e., at the end of the mixed phase after 6 fm/c), with a peak position at p^^^^ ~ 0.66 GeV, 
indeed closely resemble a thermal distribution with a temperature of about 290 MeV. Even though 
this indicates that the spectra are not fully thermalized at Tc, the change from an initial average 

{pj) = 1.66 GeV (corresponding to a "temperature" of ~920 MeV) is appreciable. The fact that 
most of rescattering occurs in the early evolution phase (the first 3 fm/c or so, cf. right panel of 
Fig. ISJ should provide favorable conditions for the build-up of elliptic flow. 

We recall that our treatment becomes unreliable toward high pr, since (i) we have neglected 
the momentum dependence of the drag and diffusion coefficients (i.e., used their values at zero 
momentum), whereas in reality they decrease with p^ (cf. Fig. [2j, and (ii) we have not accounted 
for induced gluon emission which is expected to be the main mechanism for energy loss of high- 
Pt partons within pQCD 0, 0|. Also note that transverse-flow effects, i.e., Lorentz boosts of 
the (partially) thermalized c-quarks from the comoving (thermal) frame into the Lab frame, are 
not included in Fig. El Finally, the (possibly gradual) disappearance of resonance states toward 
high temperatures has not been incorporated (whether this happens under RHIC conditions with 
initial temperatures of 2 is not clear at present; e.^., in the lattice calculations of Refs. sO, ^ 
resonance signals are still observed around these temperatures). 



IV. CONCLUSIONS AND OUTLOOK 

In the present work we have studied the role of resonant rescattering for heavy quarks in a 
Quark-Gluon Plasma at moderate temperatures. Our main assumption has been the existence of 
D-meson-like resonance states above Tc which finds support in both effective quark models and 
recent QCD lattice calculations. The underlying Lagrangian embodied both chiral and heavy- 
quark symmetry, where the latter has been essential to ensure conserved vector currents (for D* 
resonances). Bare masses and coupling constants of the model have been adjusted to render one- 
loop resummed D-meson spectral functions reminiscent to more microscopic calculations. 

Pertinent cross sections for resonant c-quark rescattering on light antiquarks in the heat bath 
have been applied within a Fokker-Planck equation to evaluate kinetic thermalization. Our main 
finding is that the introduction of resonances in the QGP leads to a substantial reduction of the 
equilibration-time scales (by a factor of ~3) as compared to estimates based on perturbative inter- 
actions. To a large extent, this difference originates from the isotropy of the angular distributions 
for resonant scattering, as opposed to mostly forward scattering in pQCD, which importantly enters 
into the (angular-weighted) transport cross section. The reduced timescales are very reminiscent to 
expected QGP lifetimes in central Au-Au collisions at RHIC. Consequently, time-evolved c-quark 
PT-spectra exhibit a marked tendency toward thermalization in the resonance picture, whereas 
they are only little affected by pQCD rescattering alone. This should have important consequences 
for the build-up of elliptic flow of charmed hadrons. 

Let us briefly discuss further ramifications and directions for future work. Clearly, the three- 
momentum dependence of the drag and diffusion coefficients, the effects of transverse flow and the 
evaluation of elliptic flow need to be addressed. Quasiparticle masses for light quarks and gluons 
should be introduced to make a closer connection to the QGP equation of state. The disappearance 
of the bound states needs to be accounted for, especially for applications at LHC with larger 
anticipated initial temperatures than at RHIC. Indeed, a key question that may eventually be 



bative contributions to, say, c + g—>c + 5 in the D-meson resonance propagator. In principle, this would lead to 
a slight renormalization of the resonance mass and width, which, however, is well inside the range of uncertainties 
of the in-medium resonance parameters. 
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answered by lattice QCD is whether I?-meson (and other) resonances in the QGP exist, and, if so, 
whether they are located above the two-quark threshold, which renders them accessible for direct 
(c+g— >c+g) scattering processes. Clearly, if this is not the case, other processes, such as c+q— >-D+q, 
need to be evaluated. The effects on secondary cc production should also be checked |42]. E.g., 
according to recent transport calculations [2^, upscaling pQCD cross sections by a factor of 3 
(to generate a significant elliptic flow) entails an increase in open-charm pairs by 40-50% over 
primordial production in central Au-Au collisions at RHIC, which is not supported by current 
PHENIX dkta 0. Resonance cross sections, due to the larger D -meson mass in the crossed 
channel, may not have this feature (or, at least, to a less extent). Finally, we recall the important 
impact that heavy-quark momentum distributions have on secondary production ("coalescence") 
of charmonium and bottomonium states. Obviously, the combined theoretical and experimental 
study of heavy-flavor probes promises a rich potential for the understanding of the complex nature 
of QGP at moderate temperatures. 
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